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ABSTRACT 

Spectral lag, which is defined as the difference in time of arrival of high- and low-energy 
photons, is a common feature in gamma-ray bursts (GRBs). Previous investigations have 
shown a correlation between this lag and the isotropic peak luminosity for long duration 
bursts. However, most of the previous investigations used lags extracted in the observer frame 
only In this work (based on a sample of 43 Swift long GRBs with known redshifts), we present 
an analysis of the lag-luminosity relation in the GRB source frame. Our analysis indicates 
a higher degree of correlation -0.82 ± 0.05 (chance probability of ~5.5 x 10“ 5 ) between 
the spectral lag and the isotropic peak luminosity, Lj S0 , with a best-fitting power-law index of 
— 1 .2 ± 0.2, such that a lag -12 . In addition, there is an anticorrelation between the source- 

frame spectral lag and the source-frame peak energy of the burst spectrum, £ P k(l + z). 

Key words: gamma-ray burst: general. 


1 INTRODUCTION 

Gamma-ray bursts (GRBs) are extremely energetic events and pro- 
duce highly diverse light curves. A number of empirical correlations 
between various properties of the light curves and GRB energetics 
have been discovered. However, the underlying physics of these 
correlations is far from being understood. 

One such correlation is the relation between isotropic peak lu- 
minosity of long bursts and their spectral lags (Norris, Marani & 
Bonnell 2000). Various authors have studied this relation using ar- 
bitrary observer-frame energy bands of various instruments (Norris 
2002; Gehrels et al. 2006; Schaefer 2007; Hakkila et al. 2008; 
Ukwatta et al. 2010c, hereafter U10). These investigations support 
the existence of the relation, however with considerable scatter in 
the extracted results. Recently, Margutti et al. (2010) investigated 
spectral lags of X-ray flares and found that X-ray flares of long 
GRBs also exhibit the lag-luminosity correlation observed in the 
prompt emission. 

*E-mail: tilan.ukwatta@gmail.com 


The spectral lag is defined as the difference in time of arrival 
of high- and low-energy photons and is considered to be positive 
when the high-energy photons arrive earlier than the low-energy 
ones. Typically, the spectral lag is extracted between two arbitrary 
energy bands in the observer frame. However, because of the red- 
shift dependence of GRBs, these two energy bands can correspond 
to a different pair of energy bands in the GRB source frame, thus 
potentially introducing an arbitrary energy dependence to the ex- 
tracted spectral lag. 

In order to explore whether the lag-luminosity relation is intrinsic 
to the GRB, it is preferable to extract spectral lags in the source 
frame as opposed to the observer frame. At least two corrections 
are needed to accomplish this: (1) correct for the time-dilation effect 
(z-correction) and (2) take into account the fact that for GRBs with 
various redshifts, observed energy bands correspond to different 
energy bands at the GRB source frame (^correction; Gehrels et al. 
2006). 

The first correction is straightforward and is achieved by multi- 
plying the extracted lag value (in the observer frame) by (1 + z) _1 . 
The second correction, on the other hand, is not so straightforward. 
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Gehrels et al. (2006) attempted to approximately correct the spec- 
tral lag by multiplying the lag value (in the observer frame) by 
(1 4- z) 033 . We note here that this correction is based on the as- 
sumption that the spectral lag is proportional to the pulse width and 
that the pulse width itself is proportional to the energy (Fenimore 
et al. 1995; Zhang et al. 2009). These approximations depend on 
clearly identifying corresponding pulses in the light curves of each 
energy land, and may be of limited validity for a large fraction 
of GREs in which the light curves are dominated by overlapping 
multipuise structures. 

Using a sample of 31 Swift GRBs, U10 found that the correlation 
coefficient improves significantly after the ^-correction is applied. 
However, this correlation does not improve further after the appli- 
cation of the ^-correction as defined by Gehrels et al. (2006). 

An alternative is to make the ^-correction by choosing two appro- 
priate energy bands fixed in the GRB source frame and projecting 
these bands into the observer frame using the relation = 

EsomceK ' 4- z). Ukwatta et al. (2010b) used this method for the first 
time to investigate the lag-luminosity relation in the source frame of 
the GRB. They selected two source-frame energy bands (100-200 
and 300-400 keV) and used background-subtracted as well as non- 
backgro md-subtracted Swift data to extract lags. Non-background- 
subtracted data were used to improve the signal-to-noise ratio for 
weak bursts. They found that the source-frame relation seems a bit 
tighter, init with a slope consistent with previous studies. Arimoto 
et al. (2010) also looked at a limited sample of High Energy Tran - 
sient Explorer (HETE)-II bursts (eight GRBs) both in the observer 
frame and the source frame, and concluded that there is no signif- 
icant effect from the redshift. However, the redshift distribution of 
their burst sample is very narrow and peaks around 1. In contrast 
to Ukwata et al. (2010b), in this study we used only background- 
subtracted data and measured the lag between source-frame energy 
bands 100-150 and 200-250 keV (the reason for selecting these 
particuliLT energy bands is described in Section 2) for a sample of 
43 Swift bursts with spectroscopic redshifts. 

In this work, we have investigated only long GRBs, i.e. bursts with 
duration greater than ~2 s. It is rather difficult to test the lag-L^ 
relation effectively for short GRBs due to a lack of spectroscopically 
measured redshifts. None of the short bursts detected so far has 
any redshift measurements obtained from a spectroscopic analysis 
of their optical afterglow. Moreover, it has been shown that short 
GRBs have either small or negligible lags (Norris & Bonnell 2006; 
Zhang e: al. 2006). According to the lag-Li*, relation, these small 
lag values imply short bursts to be highly luminous. However, based 
on the redshift measurements of their host galaxies, we can show 
that short GRBs are generally less luminous than long bursts. Hence, 
short bursts seem to not follow the lag-luminosity relation (Gehrels 
et al. 2006). 

The structure of this paper is the following. In Section 2, we 
discuss briefly our methodology for extracting spectral lags. In 
Section 3, we present our results for a sample of 43 Swift GRBs. 
We discuss our results with two candidate models in Section 4. 
Finally, in the last section (Section 5), we summarize our results 
and conclusions. Throughout this paper, the quoted uncertainties 
are at the 68 per cent confidence level. 


2 METHODOLOGY 

The Swift Burst Alert Telescope (BAT) is a highly sensitive instru- 
ment using a coded-mask aperture (Barthelmy et al. 2005). BAT 
uses the shadow pattern resulting from the coded mask to facilitate 


localization of the source. When a gamma-ray source illuminates 
the coded mask, it casts a shadow on to a position-sensitive detector. 
The shadow cast depends on the position of the gamma-ray source 
on the sky. If one knows the tile pattern in the coded mask and 
the geometry of the detector, it is possible to calculate the shadow 
patterns created by all possible points in the sky using a ray-tracing 
algorithm. Hence, by correlating the observed shadow with the pre- 
calculated shadow, one can find the location of the source. However, 
each detector can be illuminated by many sources and a given source 
can illuminate many detectors. Hence, in order to disentangle each 
sky position, special algorithms have been developed and integrated 
into the data analysis software by the Swift BAT team. 

To generate background-subtracted light curves, we used a pro- 
cess called mask weighting. The mask weighting assigns a ray- 
traced shadow value for each individual event, which then en- 
ables the user to calculate light curves or spectra. We used the 
batmaskwtevt and batbinevt tasks in ftools to generate mask- 
weighted, background-subtracted light curves, for various observer- 
frame energy bands, as shown in Table 1 . These are the energy bands 
that correspond to fixed energy bands in the source frame, i.e. 100- 
150 and 200-250 keV. These particular energy bands were selected 
so that after transforming to the observer frame they lie in the de- 
tectable energy range of the Swift BAT instrument (see Fig. 1). Even 
though the BAT can detect photons up to 350 keV, we limited the 
upper boundary to 200 keV in the observer frame. This is because 
the mask-weighted effective area of the detector falls rapidly after 
200 keV, and as a result the contribution to the light curve from 
energies greater than ^200 keV (in observer frame) is negligible 
(Sakamoto et al. 201 1 ). 

The spectral lags were extracted using the improved cross- 
correlation function (CCF) analysis method described in U10. In 
this method, the spectral lag is defined as the time delay corre- 
sponding to the global maximum of the CCF. The CCF with a delay 
index d is defined as 

CCF (d, x, y ) = ' {1) 


where x t and y L are two sets of time-sequenced data spread over N 
bins. The time delay is obtained by multiplying d by the time bin 
size of the light curves. A Gaussian curve was fitted to the CCF 
(plotted as a function of time delay) to extract the spectral lag. The 
uncertainty in the spectral lag is obtained by simulating 1000 light 
curves using the Monte Carlo technique (see U10 for more details). 

The isotropic peak luminosity (£**,) and its uncertainty for each 
GRB are obtained using the method described in U10. In essence, a 
typical GRB spectrum can be described by the Band function (Band 
et al. 1993), for the photon flux per unit photon energy using 


N(E) = 


A (ujjgv ) ° e _(2+o)£/B i ,k , E < 



( 2 ) 


which has four model parameters: the amplitude (A), the low-energy 
spectral index (a), the high-energy spectral index (ft) and the peak 
(Epic) of E?N(E) spectrum (also called the vF v spectrum, apart from 
a factor of Planck’s constant). Using these spectral parameters, the 
observed peak flux can be calculated for the source-frame eneigy 
range E\ = 1 .0 keV to £2 = 10 000 keV using 


fdos 


=L 


£i/<l+z) 

EiKl+z) 


N(E)EdE. 


( 3 ) 
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Table 1. The observer-frame energy bands and energy gaps (the energy difference between the mid-points 
of energy bands) for bursts in the sample. 


GRB 

Redshift 

Low-energv band (keV) 

High-energy band (keV) 

Energy gap (keV) 

GRB 050401 

2.899 1 

26-38 

51-64 

26 

GRB 050603 

2.821 2 

26-39 

52-65 

26 

GRB 050922C 

2.199 3 

31-47 

63-78 

32 

GRB 051111 

1.549“ 

39-59 

78-98 

39 

GRB 060206 

4.056 s 

20-30 

40-49 

20 

GRB 060210 

3.913 6 

20-31 

41-51 

21 

GRB 060418 

1.490 7 

40-60 

80-100 

40 

GRB 060904B 

0.703 s 

59-88 

117-147 

59 

GRB 060908 

1.884 9 

35-52 

69-87 

35 

GRB 060927 

5.464 1 ® 

15-23 

31-39 

16 

GRB 061007 

1.262 11 

44-66 

88-111 

45 

GRB 061021 

0.346 12 

74-111 

149-186 

75 

GRB 061121 

1.315 13 

43-65 

86-108 

43 

GRB 070306 

1.496 14 

40-60 

80-100 

40 

GRB 071010B 

0.947 15 

51-77 

103-128 

52 

GRB 071020 

2.145 16 

32-48 

64-79 

32 

GRB 0803 19B 

0.937 17 

52-77 

103-129 

52 

GRB 080319C 

1.949 18 

34-51 

68-85 

34 

GRB 080411 

1.030 19 

49-74 

99-123 

50 

GRB 080413A 

2.433 20 

29-44 

58-73 

29 

GRB 080413B 

1.101 21 

48-71 

95-119 

48 

GRB 080430 

0.767 22 

57-85 

113-141 

56 

GRB 080603B 

2.6S9 23 

27-41 

54-68 

27 

GRB 080605 

1.640 24 

38-57 

76-95 

38 

GRB 080607 

3.036 25 

25-37 

50-62 

25 

GRB 080721 

2.591 26 

28-42 

56-70 

28 

GRB 080916A 

0.689 27 

59-89 

118-148 

59 

GRB 081222 

2.770 28 

27-40 

53-66 

26 

GRB 090424 

0.544 29 

65-97 

130-162 

65 

GRB 090618 

0.540 30 

65-97 

130-162 

65 

GRB 090715B 

3.000 31 

25-38 

50-63 

25 

GRB 090812 

2.452 32 

29-43 

58-72 

29 

GRB 090926B 

1.240 33 

45-67 

89-112 

45 

GRB 091018 

0.97 1 34 

51-76 

101-127 

51 

GRB 091020 

1.710 35 

37-55 

74-92 

37 

GRB 091024 

1.091 36 

48-72 

96-120 

48 

GRB 091029 

2.752 37 

27-40 

53-67 

27 

GRB 091208B 

1.063 38 

48-73 

97-121 

49 

GRB 100621 A 

0.542 39 

65-97 

130-162 

65 

GRB 100814A 

I.440 40 

41-61 

82-102 

41 

GRB 100906 A 

1.727 41 

37-55 

73-92 

37 

GRB 1 10205 A 

2.220 42 

31-47 

62-78 

31 

GRB 11021 3 A 

1.460* 3 

41-61 

81-102 

41 


References: Batson et al. (2006); 2 Berger & Becker (2005); 3 Piranomonte et al. (2008); 4 Penprase et al. 
(2006); 5 Fynbo et al. (2009b); 6 Fynbo et al. (2009b); 7 Prochaska et al. (2006); 8 Fynbo et al. (2009b); 
9 Fynbo et al. (2009b); l0 Fynbo et al. (2009b); ll Fynbo et al. (2009b); 12 Fynbo et al. (2009b); 13 Fynbo 
et al. (2009b); 14 Jaunsen et al. (2008); l5 Cenko et al. (2007); 16 Jakobsson et al. (2007); 17 D’Elia et al. 
(2009); 18 Fynbo et al. (2009b); 19 Fynbo et al. (2009b); ^Fynbo et al. (2009b); 21 Fynbo et al. (2009b); 
^Cucchiara & Fox (2008); “Fynbo et al. (2009b); 24 Fynbo et al. (2009b); 25 Prochaska et al. (2009); 
^Fynbo et al. (2009b); ^Fynbo et al. (2009b); 28 Cuccbiara et al. (2008); 29 Chomock et al. (2009a); 
^Cenko et al. (2009); 31 Wiersema et al. (2009a); 32 de Ugarte Postigo et al. (2009); 33 Fynbo et al. (2009a); 
34 Chen et al. (2009); 35 Xu et al. (2009); 36 Cucchiara, Fox &. Tanvir (2009); 37 Chomock, Perley & Cobb 
(2009b); 38 Wiersema et al. (2009b); 39 Milvang-Jensen et al. (2010); ^O’Meara, Chen & Prochaska 
(2010); 41 Tanvir, Wiersema & Levan (2010); 42 Cenko, Hora & Bloom (2011); 43 Milne & Cenko (2011). 


The isotropic peak luminosity is defined by 


where d ' L is the luminosity distance; 


Liso — 47tfl? L /ohs» 


( 4 ) 


(i + z)c r <w 

Ho Jo y/ 4" z') 3 + 


( 5 ) 
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Figure 1. Fixed energy bands at the GRB source frame are projected to various energy bands at the observer frame, depending on the redshift. 


For the current universe, we take = 0.27, = 0.73 and the 

Hubble constant H 0 to be 70km s _l Mpc -1 (Komatsu et al. 2009). 
For more details of the L**, calculation, see U10. 

3 RESULTS 

We employed an additional 12 long bursts to the GRB sample 
(31 GRBs) that was used in U10, which increased the total sample 
to 43. This sample has redshifts ranging from 0.346 (GRB 061021) 
to 5.464 (GRB 060927), with an average redshift of ~2.0. The 
spectral information for the additional 12 bursts used in this paper 
is given in Table 2. The calculated peak isotropic luminosities, 
spanning three orders of magnitude, are given in U10 and Table 2. 

By choosing appropriate energy bands in the observer frame (ac- 
cording to the redshift of each burst), we extracted mask- weighted 
background-subtracted light curves for the selected source-frame 
energy bands 100-150 and 200-250 keV. The observer-frame en- 
ergy bands used for each burst are shown in Table 1 . Note that the 
energy gap between the mid-points of the two source-frame energy 
bands is fixed at lOOkeV, whereas in the observer frame, as ex- 
pected, this gap varies depending on the redshift of each burst (see 


Table 1). For example, in GRB 060927, this gap is 16keV and in 
GRB 061021 it is 75keV. This is in contrast to the spectral lag ex- 
tractions performed in the observer frame where this gap is treated 
as a constant. 

The extracted spectral lags for the source-frame energy bands 
100-150 and 200-250 keV are listed in Table 3. The Swift BAT 
trigger ID, the segment of the light curve used for the lag extraction 
( T -h X s and T + X E , T is the trigger time), the time binning of 
the light curve and the Gaussian curve fitting range of the CCF 
versus time delay plot (with start time and end time denoted as LS 
and LE, respectively) are also given in Table 3. Of 43 bursts in the 
sample, there are 24 bursts which have lags greater than zero. The 
remaining 19 bursts have lags either consistent with zero (16 bursts) 
or negative values (three bursts). 

For the 24 bursts which have positive lags with significance 1 o 
or greater (see Table 3), we find that the redshift-corrected lag is 
anticorrelated with Li S0 . The correlation coefficient for this rela- 
tion is —0.82 ± 0.05 with a chance probability of ^5.54 x 10“ 5 . 
The extracted correlation coefficient is significantly higher than the 
correlation coefficient (averaged over the six combinations of stan- 
dard BAT energy channels) of - — 0.68 reported in U10. Various 


Table 2. GRB redshift and spectral information. Note that uncertainties of parameters that are reported with 90 per cent confidence level have been reduced to 
ler level for consistency. 


GRB 

Peak flux 0 

V 

a 


(ergs [ ) 

Reference 

GRB 090812 

3.60 ±0.13 


n -2+0.04 

l,lw -a04 

-2J50lS;‘| 

(7.86li;|f) x 10 52 

Baumgarmer et al. (2009b); Pal* Shin et al. (2009) 

GRB 090926B 

3.20 ±0.19 

9lt} 

—0 l'l+O-O* 
u * l<5 -0.04 

—2 IS 4 " 0,3 * 
z ^°-0.31 

(5.22+q “) x 10 51 

Baumgartner et al. (2009c); Briggs (2009) 

GRB 091018 

10.30 ± 0.25 

28^° 

i < 0+0 24 

0j7 

-2.44±®!| 

(6.96+j'j|) x 10 51 

Golenetskii et al. (2009a); Markwardt et al. (2009) 

GRB 091020 

4.20 ±0.19 

47^4 

o m+o- 25 

UZU -0.25 

-1-70^1 

( 2 . 8 llg;l|) x I 0 52 

Chaplin (2009); Palmer et al. (2009) 

GRB 091024 

2.00 ±0.19 

CAA +100 

JU ^_100 

— 1 in f0 * 13 
11U -0.13 

—2 3fi +0,31 

200-0 31 

(5.56!££) x 10 s * 

Golenetskii et al. (2009b); Sakamoto et al. (2009) 

GRB 091029 

1.80 ±0.06 

6111S 

-1 46 +0,1? 

1-+O_oi7 

0 o/C+0.3] 
-2.3 o _ 031 

(1.67!°;®) x 10 52 

Barthelmy et al. (2009) 

GRB 091208B 

15.20 ±0.63 

124111 

— 1 44 +004 
1 -^—0.04 

o 20+0.29 
Z< 3 Z - 0 t .12 

( 1 . 68 lS§) x 10 52 

Baumgartner et al. (2009a);McBreen (2009) 

GRB 100621 A 

12.80 ±0.19 

951?, 

— 1 to 4 * 008 
l -' U -0.08 

-2-45ii;^ 

(2.5S!j>;£) x 10 51 

Golenetskii et al. (2010a); Ukwatta et al. (2010b) 

GRB 100314A 

2.50 ±0.13 

1061’ 


- 2 - 02 ^ 

(8.27!j;g) x 10 51 

Krimm et al. (2010); von Kienlin (2010) 

GRB 100906A 

10.10 ±0.25 

180^| 


—2 2 O +01 ^ 
- 2 . 2 U _ 013 

(4.90!i;g) x 1052 

Barthelmy et al. (2010); Golenetskii et al. (2010) 

GRB 1102Q5A 

3.60 ±0.13 

222 :* 

1 <0+0.09 
I.jZ_o 09 

—2 36* 0 * 31 
z - 3O -0.31 

oo 

^3 

X 

A 

% 

Golenetskii et al. (2011) 

GRB 110213A 

1.60 ± 0.38 

98+* 

_ i 44+0.03 

1 ^-0.03 

—2 Ifi 4 ^' 31 

Z.JO_o3i 

(3.53!‘;’’) x 10 51 

Barthelmy et al. (201 1); Foley (2011) 


a l-s peak photon flux measured in photons cm 2 s 1 in the energy range 15-150keV. 
^Peak energy, E pk , is given in keV. 
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Table 3. Source-frame spectral lag values of long duration Swift BAT GRBs. 


GRB 

Trigger ID 

T + X S (s) 

T + Xe(s) 

Bin size (ms) 

LS(s) 

LE(s) 

Lag value (ms) 

Significance 

GRB 050401 

113120 

23.03 

29.43 

64 

-2.00 

2.00 

310 ± 145 

2.14 

GRB 050603 

131560 

-3.83 

3.08 

16 

-0.40 

0.40 

-16 ±21 

-0.76 

GRB 050922C 

156467 

-2.70 

2.94 

16 

-1.00 

1.00 

136 ±68 

2.00 

GRB 051111 

163438 

-6.96 

28.62 

64 

-4.00 

4.00 

333 ±251 

1.33 

GRB 060206 

180455 

-1.29 

8.18 

16 

-2.00 

2.00 

86 ± 111 

0.77 

GRB 060210 

180977 

-3.37 

5.08 

128 

-4.00 

4.00 

658 ±259 

2.54 

GRB 060418 

205851 

-7.66 

33.04 

64 

-2.00 

2.00 

-110 ± 106 

— 1.04 

GRB 060904B 

228006 

-1.97 

10.32 

512 

—6.00 

6.00 

124 ±436 

0.28 

GRB 060908 

228581 

-10.91 

3.68 

32 

-2.00 

2.00 

78 ± 124 

0.63 

GRB 060927 

231362 

-1.69 

8.04 

32 

-1.00 

1.00 

18 ±75 

0.24 

GRB 061007 

232683 

23.86 

65.08 

4 

-0.20 

0.20 

52 ±22 

2.36 

GRB 061021 

234905 

-0.46 

14.64 

512 

-4.00 

4.00 

-430 ± 975 

-0.44 

GRB 061121 

239899 

60.44 

80.66 

4 

-0.20 

0.20 

22 ± 10 

2.20 

GRB 070306 

263361 

90.00 

118.42 

32 

-4.00 

2.00 

-362 ± 247 

-1.47 

GRB 071010B 

293795 

-1.70 

17.24 

64 

-2.00 

2.00 

404 ± 159 

2.54 

GRB 071020 

294835 

-3.22 

1.14 

4 

-0.20 

0.40 

35 ±13 

2.69 

GRB 080319B 

306757 

-2.85 

5151 

4 

-0.10 

0.14 

23 ±6 

3.83 

GRB 0803 19C 

306778 

-0.77 

13.31 

32 

-1.00 

1.00 

174 ±91 

1.91 

GRB 080411 

309010 

38.46 

48.45 

4 

-0.50 

0.50 

116 ±25 

4.64 

GRB 080413A 

309096 

-0.42 

9.05 

8 

-1.00 

1.00 

107 ± 59 

1.81 

GRB 0804 13B 

309111 

-1.44 

4.96 

32 

-1.00 

1.00 

115 ±50 

2.30 

GRB 080430 

310613 

-1.24 

12.84 

256 

-4.00 

4.00 

91 ±431 

0.21 

GRB 080603B 

313087 

-0.54 

5.10 

16 

-1.00 

1.00 

5 ±59 

0.08 

GRB 080605 

313299 

-5.46 

15.53 

8 

-0.20 

0.20 

35 ±18 

1.94 

GRB 080607 

313417 

-6.13 

12.05 

8 

-0.50 

0.50 

26 ±30 

0.87 

GRB 080721 

317508 

-3.39 

8.64 

64 

-2.00 

2.00 

-86 ± 110 

-0.78 

GRB 080916A 

324895 

-2.66 

39.58 

128 

-2.00 

4.00 

585 ± 214 

2.73 

GRB 081222 

337914 

-0.80 

15.58 

4 

-1.00 

1.00 

227 ± 51 

4.45 

GRB 090424 

350311 

-0.94 

4.95 

16 

-0.20 

0.20 

14 ±14 

1.00 

GRB 090618 

355083 

46.01 

135.35 

8 

-2.00 

2.00 

267 ± 72 

3.71 

GRB 090715B 

357512 

-4.80 

21.06 

16 

-2.00 

3.00 

275 ± 155 

1.77 

GRB 090812 

359711 

-6.93 

41.20 

256 

-6.00 

6.00 

-22 ±202 

-0.11 

GRB 090926B 

370791 

-22.00 

36.00 

512 

-10.00 

8.00 

746 ±627 

1.19 

GRB 091018 

373172 

-0.28 

2.92 

64 

-2.00 

1.00 

143 ± 297 

0.48 

GRB 091020 

373458 

-2.54 

13.84 

128 

-3.00 

2.00 

-187 ± 177 

-1.06 

GRB 091024 

373674 

-9.58 

27.29 

512 

-10.00 

10.00 

912 ±604 

1.51 

GRB 091029 

374210 

-4.03 

38.98 

256 

-10.00 

10.00 

-112 ±395 

-0.28 

GRB 091208B 

378559 

7.66 

10.61 

64 

-1.00 

1.00 

105 ± 66 

1.59 

GRB 100621A 

425151 

-6.79 

40.31 

256 

-3.00 

3.00 

1199 ±311 

3.86 

GRB 1008 14A 

431605 

-4.40 

29.39 

256 

-4.00 

4.00 

862 ± 147 

5.86 

GRB 100906 A 

433509 

-1.49 

26.16 

128 

-2.00 

2.00 

105 ±79 

1.33 

GRB 11 0205 A 

444643 

118.89 

293.99 

64 

-1.00 

1.00 

-29 ± 52 

-0.56 

GRB 110213A 

445414 

-3.42 

5.29 

512 

-3.00 

3.50 

602 ± 746 

0.81 


Table 4. Correlation coefficients of the lag-luminosity 
relation. 


Coefficient type 

Correlation coefficient 

Null probability 

Pearson’s r 

-0.82 ± 0.05 

5.54 x 1<T 5 

Spearman’s r s 

-0.70 ± 0.06 

1.49 x 10- 4 

Kendall’s r 

-0.50 ± 0.05 

6.63 x 10 -4 


correlation coefficients of the relation are shown in Table 4, where 
uncertainties in the correlation coefficients were obtained through 
a Monte Carlo simulation utilizing uncertainties in and the 
lag values. The null probability that the correlation occurs due to 
random chance is also given for each coefficient type. 

Fig. 2 shows a log-log plot of isotropic peak luminosity versus 
redshift-corrected spectral lag. The solid line shows the following 


best-fitting power-law curve: 

log ( = (54.7 ± 0.4) - (1.2 ± 0.2) log (6) 

yerg s _1 y 1+z 


Since there is considerable scatter, the uncertai nties of the fit pa- 
rameters are multiplied by a factor of \/x 2 /ndf = .36/22 = 

1.96. The dash lines indicate the estimated lcr confidence level, 
which is obtained from the cumulative fraction of the residual dis- 
tribution taken from 16 to 84 per cent. 

The best-fitting power-law index (— 1 .2 ± 0.2) is consistent with 
observer-frame results obtained by Norris et al. (2000) (~— 1.14) 
and the average power-law index of — 1 A ± 0.3 reported in U10. 
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Figure 2. The spectral lags between the source-frame energy range bands 100-150 and 200-250 keV and the isotropic peak luminosity are plotted in a log-log 
plot. 


4 DISCUSSION 

4.1 Spectral lags: observer frame versus source frame 

U10 extracted spectral lags in fixed energy bands in the observer 
frame, and in this work for the same sample of 3 1 bursts we extracted 
lags in fixed energy bands in the source frame. In the observer- 
frame case, there are four energy channels [canonical BAT energy 
bands: channel 1 (15-25 keV), 2 (25-50 keV), 3 (50-100 keV) and 
4 ( 1 00-200 keV)], thus six lag extractions per burst. It is interesting 
to study to what degree these different lags correlate with source- 
frame lags (between fixed source-frame energy channels 100-150 
and 200-250 keV). In Fig. 3, we show all combinations of observer- 
frame lags as a function of source-frame lags. The red data points 
show lags with the time-dilation correction due to cosmological 
redshift, and black data points show lags without the time-dilation 
correction. From Fig. 3, it is clear that all plots show some corre- 
lation both in the time-dilation-corrected (shown in red) and time- 
dilation-uncorrected (shown in black) cases. We note that the corre- 
lation coefficients are greater than 0.5 in time-dilation-uncorrected 
cases where BAT channel 1 is involved in the lag extraction. In the 
time-dilation-corrected case, all plots show correlation coefficients 
greater than 0.5 except for the lag 43 plot. Despite these moderate 
correlation coefficients, the large scatter seen in these plots indi- 
cates that the observer-frame lag does not directly represent the 
source-frame lag. 


4.2 Lag-Lfco relation: observer frame versus source frame 

There are two important changes in the lag-luminosity relation 
which may occur when going from fixed observer-frame energy 
bands to fixed source-frame energy bands: a change in the power- 
law index, and a change in the dispersion of the data measured by the 
correlation coefficient. Table 5 summarizes these two parameters for 

© 2011 The Authors, MNRAS 419, 614-4523 
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various energy bands both in the observer frame and in the source 
frame. 

In the observer frame, the power-law index varies from ~0.6 to 
^ 1 .8, with mean around 1 .3. In the source frame, the index changes 
from 0.9 to 1.23 with a mean of ~ 1.1. Meanwhile, the correlation 
coefficient varies from 0.60 to 0.79 in the observer frame, and in 
the source frame it changes from 0.76 to 0,90. Hence, according to 
Table 5, the source-frame lag-L^ relation seems to be tighter than 
the observer-frame case with a slope closer to 1. 

4.3 Spectral lag-f pk relation 

Now we investigate the relation between source-frame spectral lag 
and source-frame average peak energy [E^(\ + z)] of the burst 
spectrum. In Fig. 4, we plotted £ pk (l + z) as a function of source- 
frame lags. There is a correlation between these two parameters 
with a correlation coefficient of —0.57 ± 0.14. Various correlation 
coefficients of the relation are shown in Table 6, with uncertainties 
and null probabilities. 

The best fit is shown as a dashed line in Fig. 4, yielding the 
following relation between £ pk (l + z) and lag/(l + z): 

108 ( Ep tev ~~) = (3J ±01)_ (0 ' 56 ± 006)108 (7) 

The uncertainties in the fitted parameters are expressed with the 
factor of x/xVndf = y/Wl\/22 * 1.18. 

According to equation (6), oc [lag/(l 4- z)] -12 . From the 
Yonetoku relation, we know that Lj S oOc[£ pk (l + z)] 20 (Yone- 
toku et al. 2004). Hence, from these two relations, we expect 
to see a correlation between £ pk (l + z) and lag/(l + z) such as 
£pk(l+z)oc [lag/(l +z)]-°‘ 6 . 

The best-fitting slope of 0.56 ± 0.06 is consistent with the ex- 
pected slope of ~0.6 based on the source-frame lag-luminosity and 
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Figure 3. All combinations of fixed observer-frame energy channel [canonical BAT energy bands: channel 1 (15-25 keV), 2 (25-50 keV), 3 (50-100 keV) and 
4 (100-203 keV)] spectral lag values as a function of fixed source-frame energy channel (between 100-150 and 200-250 keV) lag values. Black and red data 
points and labels correspond to redshlft-uncorrected and redshift-corrected cases, respectively. The blue dashed line corresponds to the equality line of the two 
parameters in each panel. 


the Yonetoku relation. However, note that the correlation coefficient 
is significantly smaller than the coefficient for the lag-luminosity re- 
lation. This lower degree of correlation may be suggestive of bright- 
ness and detector-related selection effects that have been noted in 
the literature (Butler et al. 2007) for the Yonetoku relation. 

4.4 Some models for spectral lags 

U10 and this work have provided more evidence for the existence 
of the lag-luminosity relation based on a sample of Swift BAT 


GRBs with measured spectroscopic redshifts. This analysis calls 
for a physical interpretation for spectral lag and a lag-luminosity 
relation. In the literature, several possible interpretations have been 
discussed (Dermer 1998; Salmonson 2000; Ioka & Nakamura 2001; 
Kocevski & Liang 2003; Qin et al. 2004; Schaefer 2004; Ryde 2005; 
Shen, Song & Li 2005; Lu et al. 2006; Peng et al. 2011). 

One proposed explanation for the observed spectral lag is the 
spectral evolution during the prompt phase of the GRB (Dermer 
1998; Kocevski & Liang 2003; Ryde 2005). Due to cooling effects, 
£ pk moves to a lower energy channel after some characteristic time. 


Table 5. Observer-frame and source-frame slopes and correlation coefficients of the lag-Li s0 relation. Conservative 10 per cent 
uncertainty is assumed for cases without uncertainties. 


Energy bands 

Frame 

Slope 

Couelation coefficient 

Number of GRBs 

Reference 

(0.3-1), (3-10) keV 

Observer 

0.95 ± 0.23 

_ 

9 

Margutti et al. (2010) 

(6-25), (50-400) keV 

Observer 

1.16 ±0.07 

—0 7Q+° 16 
U>/V -0.05 

8 

Arimoto et al. (2010) 

(15-25), (25-50) keV 

Observer 

1.4 ±0.1 

-0.63 ± 0.06 

21 

U10 

(15-25), (50-100) keV 

Observer 

1.5 ±0.1 

-0.60 ± 0.06 

28 

U10 

(15-25), (100-200) keV 

Observer 

1.8 ±0.1 

-0.67 ± 0.07 

27 

U10 

(25-50), (56 100) keV 

Observer 

1.2 ±0.1 

-0.66 ± 0.07 

27 

U10 

(25-50), (100-200) keV 

Observer 

1.4 ±0.1 

-0.75 ±0.07 

25 

U10 

(25-50), (100-300) keV 

Observer 

1.14 ±0.1 

- 

6 

Norris et al. (2000) 

(25-50), (100-300) keV 

Observer 

0.62 ±0.04 

-0.72 ±0.07 

6 

Hakkila et al. (2008) 

(50-100), (100-200)keV 

Observer 

1.4 ±0.1 

-0.77 ± 0.08 

22 

U10 

(20-100), (100-500) keV 

Source 

1.23 ±0.07 

-0.90*$ 

8 

Arimoto et al. (2010) 

(100-200), (300-400) keV 

Source 

0.9 ± 0.1 

-0.76 ±0.06 

22 

Ukwatta et al. (2010b) 

(100-150), (200-250) keV 

Source 

1.2 ±0.2 

-0.82 ± 0.05 

24 

This work 
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Lag/(1+z) msec 

Figure 4. The source-frame peak energy [£ pk (l 4- z)] versus source-frame spectral lags. The energy bands, 100-150 and 200-250 keV, corresponding to the 
lag extractions are shown in hashed red bands on the plot. 


Table 6. Correlation coefficients of the lag-£ pk relation. 


Coefficient type 

Correlation coefficient 

Null probability 

Pearson’s r 

—0.57 ±0.14 

4.83 x 10“ 3 

Spearman’s r 5 

-0.50 ±0.12 

1.36 x 10~ 2 

Kendall’s r 

-0.37 ±0.14 

1.18 x 10“ 2 


When the peak energy (E^) moves from a higher energy band to a 
lower energy band, the temporal peak of the light curve also moves 
from a higher energy band to a lower one, which results in the 
observed spectral lag. In a recent study, Peng et al. (2011) sug- 
gest that spectral evolution can be invoked to explain both positive 
and negative spectral lags. Hard-to-soft evolution of the spectrum 
produces positive spectral lags, while soft-to-hard evolution would 
lead to negative lags. In addition, these authors also suggest that 
soft-to-hard-to-soft evolution may produce negative lags. 

A schematic diagram showing a hard-to-soft scenario is depicted 
in Fig. 5. Initially, Epk of the spectrum is in the high-energy band, 
which results in a pulse in the light curve of the high-energy band. 
Then £ pk moves to the lower energy band resulting in a pulse in the 
low-energy light curve. The temporal difference between the two 
pulses in the light curves would then be a measure of the cooling 
time-scale of the spectrum. 

If this were the only process that caused the lag, then in a simple 
picture one would expect the source-frame average £ pk to lie within 
the two energy bands in question. According to Fig. 4, for the 
majority of bursts the source-frame E pk lies outside the energy band 
100-250 keV, indicating that the simple spectral evolution scenario 
described above may not be the dominant process responsible for the 
observed lags. However, it is worth noting that a pulse in a specific 
energy band may not always mean that the E pk is also within that 
energy band. There are other issues associated with this model: (1) 
the calculated cooling times based on simple synchrotron models 
are, in general, relatively small compared to the extracted lags and 
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Figure 5. The time evolution of the £ pk across energy bands may cause the 
observed spectral lags in GRBs. 


(2) short bursts which exhibit considerable spectral evolution do not 
show significant lags. 

Another model that purports to explain spectral lags is based on 
the curvature effect, i.e. a kinematics effect due to the observer 
looking at increasingly off-axis annulus areas relative to the line 
of sight (Salmonson 2000; Ioka & Nakamura 2001; Dermer 2004; 
Shen et al. 2005; Lu et al. 2006). Fig. 6 illustrates how the spectral 
lag could arise due to the curvature effect of the shocked shell. Due 
to a smaller Doppler factor and a path difference, the radiation from 
shell areas which are further off-axis will be softer and therefore 
lead to a lag. As with spectral evolution models, there are difficulties 
associated with the curvature models too. These kinematic models 
generally predict only positive lags. As can be seen from Table 3, 
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Figure 6. Spectral lags could arise due to the curvature effect of the shocked shell. At the source, the relativistically expanding shell emits identical pulses from 
all latitudes. However, when the photons reach the detector, on-axis photons get boosted to higher energy (hard). Meanwhile, off-axis photons get relatively 
smaller boost and travel longer to reach the detector. Thus, these photons are softer and arrive later than the on-axis photons. 


some of the measured lags are negative, and therefore these lags 
present a real challenge for the simple curvature models. 

It is possible that spectral lags are caused by multiple mecha- 
nisms. Peng et al. (2011) investigated spectral lags caused by in- 
trinsic spectral evolution and the curvature effect combined. They 
showed that the curvature effect always tends to increase the ob- 
served spectral lag in the positive direction. Even for cases with 
sofit-to-hard spectral evolution, when the curvature effect is intro- 
duced lags become positive. Hence, they predict that the majority 
of measured spectral lags should be positive, which is consistent 
with the findings of this work and U10. 


5 SUMMARY AND CONCLUSION 

We have investigated the spectral lag between 100-150 and 200- 
250 keV energy bands at the GRB source frame by projecting these 
bands to the observer frame. This is a step forward in the investiga- 
tion of lag-luminosity relations since most of the previous investi- 
gations used arbitrary observer-frame energy bands. 

Our analysis has produced an improved correlation between spec- 
tral lag (t) and isotropic luminosity over those previously reported 
with the following relation: 

1 +z 

We also find a modest correlation between the source-frame spec- 
tral lag and the peak energy of the burst, which is given by the 
relation 

r/ms [“ £ P k(l 4- £)/(keV) 

14 z [ 10 37 

Finally; we mentioned two simple models and noted their limita- 
tions in explaining the observed spectral lags. 



W(ergs~') 

10 547 
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